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In a process analogous to photon-photon creation via the dynamical Casimir effect (DCE), after 
separating Maxwell's equations using Hertz vectors in a generalized Lorenz gauge, and which allows 
for both electric and magnetic time-dependent media, we discuss how surface plasmon polaritons 
(SPPs) can be created out of vacuum via the time-dependent variation of the medium. We find 
that vacuum excited SPPs dominate the photon creation rate for TM modes in a dielectric slab, 
modulated in time inside a cylindrical cavity The possible experimental detection of vacuum SPPs 
is also discussed. 
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Particle creation via the Schwinger-effect in static 
electric fields [1], in expanding universes [2] or from 
black hole evaporation [3] all have yet to be confirmed. 
However there is a related effect: the dynamical 
Casimir effect (DCE), first discussed by Moore [4] who 
showed that pairs of photons can be created in a cavity 
if one of the walls moves non-adiabatically. The num- 
ber of photons created during a given number of para- 
metric oscillations is proportional to sinh^ {2ujt v / c) ^ 
e.g., see [5], where v is the wall velocity and c is the 
speed of light. To overcome the fact that the mechan- 
ical properties of the material usually imply v/c <C 1, 
there have been proposals other than mechanical os- 
cillations. Calculations show that modulating a di- 
electric medium using a laser leads to particle cre- 
ation by varying the optical path length of the cav- 
ity, e.g., see [6]. There are also related experiments 
underway in three-dimensions using centimeter-sized 
(microwave) cavities [7], where a laser is instead used 
to modulate the surface conductivity. Other meth- 
ods use illuminated superconducting boundaries [8] 
and recently time varied inductance effects in one- 
dimensional quantum circuits have shown vacuum 
squeezing [9, 10]. Rotating analogs have also been 
investigated [11]. 

In this letter we discuss the possibility of the cre- 
ation of vacuum excited surface plasmon polaritons 
(SPPs). In Fig. 1 the general idea is sketched, where 
a pulsed laser of an appropriate frequency can be used 
to vary the time dependence of the dielectric. The 
cavity is not per se required; however, a low tem- 
perature cavity will suppress thermal excitations and 
allows for comparison with calculations for photon- 
photon rates [6]. The fact that SPPs can be excited 
from vacuum fluctuations besides photons is much like 
the role SPPs play in the static Casimir force [12] 
with a metal/insulator/metal (MIM) heterostructure. 
We essentially generalize this to the dynamical case 
at first to the simpler case of a single metal/insulator 
(MI) interface and show that the time modulation of a 
dielectric leads not only to two-photon processes, but 
to vacuum excited SPPs that are comparable to the 
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FIG. 1. A pulsed laser uniformly irradiates (via a lens) a 
dielectric slab ( region 1) of thickness a and radius R in a 
superconducting cylindrical cavity (dotted lines) of length 
L. Vacuum excited SPPs (depicted by - + - + - +) can be 
detected by a photon scanning tunneling microscope placed 
in region 2. Alternatively a phase matched prism can be 
attached above the interface (in region 2) to detect surface 
plasmon resonances. 



photon creation rate. 

Maxwell's equations can be derived from the follow- 
ing Lagrangian: 
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mf-::m\t)<^' (1) 



(so = Mo = !)• In the above we assume that the elec- 
tric permittivity and magnetic permeability are time 
dependent, but constant in space: V/i = Ve = 0. ^ 
represents a TM field with generalized Neumann BCs 
in a cavity and the TE case (swapping £ ^ /i) is rep- 
resented by ^ with Dirichlet BCs, e.g., see [13]. This 
Lagrangian is useful because the standard canonical 
Hamiltonian can be constructed [14], where the mass 
term, {t) , represents the coupling of light to a time- 
dependent boundary (set to = in what follows). 

A convenient way to separate Maxwell's equations 
is using Hertz vectors; developed by Nisbet [15] for 
non-dispersive inhomogeneous media. However, here 
we generalize to the case of a constant isotropic, but 
time-dependent medium. It is then possible to show 
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that Maxwell's equations separate as 

sit)dtitJ.it)dtU,) - v^He = , 
fi{t)dt{e{t)dtnm) - v^n™ = , 



(2) 



where we use a generalized Lorenz gauge: 
IJ.{t)dt{s{t)Ao) + V • A = 0, cf. [15, 16]. In the 
above we have assumed both a zero permanent 
polarization and magnetization (Pq = Mq = 0) as 
well as zero bulk charges and currents (p = 0, J = 0), 
although these can also be included in the Hertz 
method. Note the Lagrangian in Eq. (1) leads to 
the equations of motion for Ilg in Eq. (2) (11^ is 
obtained by swapping ja ^ e in Eq. (1)). This 
approach generalizes other work [6, 16] which con- 
sidered only time-dependent £ or /i. Further work 
for non-dispersive, inhomogenous, conducting and 
time-dependent media: e(r, t) and /i(r, t), will be 
presented elsewhere. 

Before discussing quantized SPPs we first need to 
find the classical solutions for a single interface (be- 
tween two media) that lead to SPPs and are usually 
only discussed for TM modes [17]. With that goal 
in mind, we write the electric and magnetic fields in 
terms of the Hertz vectors: 

E = 1 V X (V X He) - /ioV X dtUm , 

dU 

B = /iV X ^ + /io V X (V X n^) , (3) 

allowing one to easily isolate TE and TM modes. 
For example, TM modes are defined by the parts 
Etm,Btm coming from Xlg with z • B = 0, where 
a convenient choice of Hertz vectors are: Xlg = 
$ z , = ^ z. For TM modes we obtain 

Etm = -didz^ei + -d2dz^e2 , 

e £ 



(4) 



with a similar expression for TE modes. These gener- 
alize the time-independent cases found, e.g., in [18]. 

In what follows take two half spaces in the z- 
direction, where region 1 (the semiconductor slab) is 
to be metallic (M) during irradiation and region 2 an 
insulator (I) like air, creating an MI interface (later 
on we will discuss layered IMI and MIM heterostruc- 
tures). In our proposed set up £i{t) varies from ± 
values and £2 remains constant (although for now it 
will be left more general). To make explicit the utility 
of the Hertz vector method we shall consider the ra- 
dial propagation of SPPs in a cylindrical cavity with 
coordinates (p, 6>, z), sectional radius p = R and length 
L, see Fig. 1. 

To allow for time-dependent mode functions in 
either region we assume an 'instantaneous' set of 
mode functions [19]: $(r,t) = Qm(^)^m(r; t) 
with m = (n,p, m), where t becomes a parame- 
ter: (p{r,t) (p{r;t) and dz e{t)(prn{T^;t)(pri{T^;t) = 
(^m,^n) = ^mn (orthouormal) and satisfy the wave 
equation: VVm(r; t) + £(t)/i(t)a;^(t)(/?m(r; t) = 0. 
If we substitute this into the Lagrangian in Eq. 



(1) using orthonormality of the mode functions and 
defining the conjugate momentum as n(r, t) = 
^(^) Xlm ^m(^)v^m(r; t) via a Legendre transform we 
then find an 'effective' Hamiltonian [20]: i^eff = 

Em {Pl+^lmi,) + EmnPmQnMmnit), wheVB 

Pm = Qm — -MmnQn ^ud the iutcrmodc coupling ma- 
trix is A^mn = Jq dz e{t)(pm{^'',t)dt(pn{^'',t)' Howcvcr, 
for SPPs, using the definition in Eq. (5) below, it 
is easy to show that there are no intermode coupling 
terms: A^mn = 0, like for a uniform dielectric filling 
the entire cavity (a = L) [19]. Promoting these to 
operators, the ETCRs then become [Qm, Ai] = ^^^mn, 
which are then equivalent to [4>m,IIn] = ^^mn when 
we quantize each degree of freedom (<l>, ^). 

To investigate SPPs for a single interface the ansatz: 



sp/ i\ _ J ^np 



(x±), 



z < a; £1, /ii , 
z > a; £2, /i2 , 
(5) 

leads to the following dispersion relations: 
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where Eq. (2) was used in each region. For SPPs, the 
transverse modes are essentially unbounded and are 
simply given by rnp(x^) = e*^^^, where s = (p, 6>). 
In cylindrical coordinates the transverse Laplacian is 
defined by V^rk^ = — k^rk^ with eigenvalue k^. 
For the photon branch the modes are bounded where 

rnp(xx) = i?j^/i(a^„p) Jn (^np^)e''^^ and Xnp IS the 
pth root of Jn(^) = [21]. Note these mode functions 
are orthonormal: (r^n^^np) = 

Standard boundary conditions at an interface: 
(D2 - Di) • z = and z X (E2 - El) = [21] then im- 
ply Ax = A2 and i^i{t)/£i{t) + i^2{t)/£2{t) = which 
requires two dielectrics to be of opposite sign to gen- 
erate SPPs [17]. Eliminating the z-dependent hCi we 
find the following 'electric' dispersion relation: 



£l£2 



£lP2 - £21^1 



(J) 



£1 ^ £2 \ ^1 ~ ^2 

Note the symmetry in reversing regions 1 and 2. 
Then with pi = 112 we get the standard result 

{^yKY — ^^^(^ + ^2/^1), where here we allow for 
time-dependent dielectrics. 

It is also possible to show that magnetic SPPs exist 
for TE modes. Using the TE components of the Hertz 
vectors and using an 'instantaneous' equation like (5) 
for V^m(r;t) along with (B2 — Bi) • z = and z x 
(H2 — Hi) = lead again to Ai = A2 but now with 
Ki{t)/iii{t) + n2{t)/ii2{t) = 0. Hence SPPs do exist 
for TE modes as long as for example, /ii < 0, /i2 > 0, 
which can be achieved using split ring resonators, e.g., 
see [17]. Finally using 11^ in Eq. (2) leads to the 
'magnetic' dispersion relation: 



P1P2 



Pl£2 - P2S1 



(8) 



Ml + M2 V Ml - M2 

This result can be obtained from the 'electric' sector 
by swapping £i ^ jii and when £1 = £2 the result 
simplifies to (c^^)^ = ^^(1 + M2/M1). 
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To compare vacuum excited SPPs with some exper- 
imental proposals for photon creation using semicon- 
ductor slabs, e.g., see [5], we will consider TM modes. 
These have the following orthonormal mode functions 
in terms of the photon Hertz scalar: 



^icos {kimz)r np{x_i), 0< z <a, 
A2 cos (k2m{L - z))rnp{y.i_), a< z < L. 

(9) 

where (using the same TM interface conditions as be- 
fore) we find following transcendental equation 



/ciTO tan(A;iTOa) k2m taT^ihmla - L]) 



S2{t) 



(10) 



for the eigenvalues. This agrees with the result in [6], 
but can be derived with the minimum of effort using 
Hertz vectors and generalized to arbitrary transverse 
section. Note that the photon dispersion relation (in 
this case for a for a cylindrical section) at any given 
time in regions i = 1,2: 
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must be equal at the interface implying that ^(A^i 

{^f) = ^(^2m+ (^)^)- Both this constraint and 
the eigenvalue relation, Eq. (10), must be simultane- 
ously satisfied [6]. 

The number of particles created can be obtained 
using the Heisenberg formalism, where the quan- 
tum field operator expansion is defined in terms 
of annihilation and creation operators as 4>(r, t) = 
Em [^m^m(r,^) + aiS'^^iS*(r,^)] • Bcforc laser irradi- 
ation (t < 0) we have a stationary solution u^^{r^t) = 

^7^=^m(r; to), while for < t < ti (during irra- 



diation) we again expand in an 'instantaneous' ba- 
sis: u'^{r,t) = EmQm(t)v^m(r;t) with initial con- 
ditions Qm(0) = -7== and Qm(0) de- 
rived from continuity of and at t = 0. When 
this is substituted into Eq. (2) and using orthonor- 
mality of (pm{'^'',t) we find an equation for an un- 
coupled set of time-dependent harmonic oscillators: 



0. 



For times 
functions in 



mn — 

t > ti we expand the mode 
terms of a new stationary eigen- 



out 



frequency cj^: u 
where the two sets of mode 
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functions 



^m(r;ti 



lated by a Bogolyubov transformation: 



t(r,t) 



Em [«nm^iS(r,t) +/5nm^iS*(r,t)] , whcrC 

((^mS<)) and f3^^ = -((^^S^S [<]*)) and the 
scalar invariant product is defined ((0,?/^)) = 

~^ /cavity ^^^(^^* ~ ^^*)- -^mn = and as- 
suming the background field (the laser) leads to 
shifts in frequency near to parametric resonance then 

^m(^) ^ ^Om + Acjgm = ^Om(l + ^ sin(2cJ0mt) ) , 

where (x;om is the driving frequency. In the late time 
limit we then have enhancement A/'m = \Pmm{t)\'^ = 
sinh^ {ujOmi^t/2) [22]. 

For simplicity consider /ii = /i2, then the time- 
dependent 'electric' SPP dispersion relation, see below 



Eq. (7), for £2 a constant and £min (-max 
with Smin < gives 1 + I ^ ^ - sin(cjomt), 

where x is an overall time-independent frequency 
shift. X = ~1 then generates particles with enhance- 
ment factor 



(12) 



which is one of the main results of this letter. We 
can now compare this to the cjqii TM mode (the low- 
est frequency cylindrical mode [13]) where in the limit 
of a <C L the photon eigenvalues shift by Aa;^(t) = 

2x'^ 

'"'Ji I _L J 

L2 



£2 



1 



o 



. Assuming the coupling 



i?2£2 L 

matrix A^mn = which can be achieved for certain 
cavity geometries [6], then the parametric enhance- 
ment for the photon branch is achieved by choosing 
— 1 ~ X ~ 1 + sin(cjom^) and leads to (for 

X = 1): = sinh^ . For a cylindrical 

cavity the SPP creation rate dominates the photon 
rate if k'j_ ^ {x'^p/ R'^){a/L)^ which for example with 
a radius R = 2.5 cm and length L = 10 cm, then if 
a/L - C>[10-'^] and xqi = 2.4048 implies kl > 1/25 
or k± ^ 1/5 which is easily achieved for SPPs of mi- 
cro to nanometer wavelength, cf. Agp = 27r/Re[/cx] 
(we discuss the propagation length L = (2Im[/c^])~^ 
later). 

To realistically modulate the permittivity, £i(t), we 
could use an appropriately doped semiconductor (with 
two well defined energy levels within the band gap) via 
Rabi oscillations, e.g, see [23]. However, in this case, 
Smin > (implying frequency shifts x > 0) and there 
are no SPP solutions. For laser pulses with an energy 
{hX/c) above the band gap, a time varying bulk con- 
ductivity, p(t), would be generated leading to a modu- 
lated permittivity with ^min < (shifts x < 0). Such a 
conductivity can be included in Maxwell's equations, 
Eq. (2), using J = (j(t)E and hence this only affects 
TM modes (surface charges, a{t) are present for both 
TE and TM modes [13]). An alternative using only 
positive definite modulations for £i{t) using IMI het- 
erostructures, is discussed later. 

Possible ways to detect vacuum excited SPPs would 
be to use near field microscopy with a photon scan- 
ning tunneling microscope, e.g., see [17], where the 
microscope is placed on the opposite vacuum/air side 
of the dielectric slab, see Fig. 1. Usually, to gener- 
ate SPPs a monochromatic light source is sent into a 
prism placed above an interface (region 2) with total 
internal reflection at angle where their presence is 
then verified by finding a decrease in emitted power 
at Oi. However, the time reversal of this process is 
equivalent to the creation of vacuum excited SPPs and 
hence might lead to the intriguing result that a signal 
can be detected from an empty prism during the time 
modulation of Si. We should of course require that the 
pulsed laser itself does not generate SPPs, as can be 
arranged by uniformly irradiating the dielectric slab 
at a right angle of incidence, see Fig. 1. 

It is also possible to consider more complicated 
multilayer regions such as insulator/metal/insulator 
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(IMI) heterostructures or MIMs (made of three re- 
gions 1,2, and 3, respectively). It is easy to verify 
(when 111 = 112) that our results follow through as in 
the standard case, but now £1, for example, is time 
dependent: e'^^-^ = /^2/.2+/.i/.i(t) ^ ^ 

particular, IMIs develop long range SPPs for 'even' 
mode solutions [17] when the metal thickness a ^ 
and may allow for easier detection of vacuum excited 
SPPs. This is also shows that SPP solutions are possi- 
ble with a modulation of ei positive definite, because 
region 2 (metal) has S2 < 0, as we mentioned earlier. 

We also discussed how SPPs can occur for TE 
modes if we assume: /i^ < (a metamaterial with 
/ii < would simultaneously allow for both TE and 
TM mode SPPs). Our results for parametric enhance- 

ment go through the same: = sinh — , 

^ ^ my uJOm 2/i2 J ' 

with /i2 replacing 62- We expect this rate to also be 
comparable to the TE photon-photon rate if fii were 
varied somehow in time at GHZ frequencies. It might 
be of further interest to try to design experiments in 
centimeter sized cavities using split ring resonators 
and wire rods that are then modulated in time. This 
may allow for easier detection by precisely controlling 
the SPP wavelength, Agp. 

In summary, in this letter we discussed how SPPs 
might be excited from vacuum for the case of a di- 
electric slab which behaves as an insulator /insulator 
(II) interface and changes to an MI interface during 
irradiation. We separated Maxwell's equations us- 
ing the Hertz vectors approach for time-dependent 
media finding both electric and magnetic SPP solu- 
tions. For parametric oscillations of the dielectric in- 
side a superconducting microwave cylindrical cavity 
(cooled to suppress thermal excitations) our calcula- 
tions show that vacuum excited SPPs can dominate 
over the photon-photon rate. The results for the pho- 
ton creation rate in a cylindrical cavity generalizes 
those found for the rectangular case [6]. 

We used an 'instantaneous' basis, which treats time 
as a parameter and might warrant some concern re- 
garding quantization. Firstly, we are not consider- 
ing moving boundaries which lead to non-separable 
Hilbert spaces; however, these can be remedied by 
a coordinate transformation to a stationary prob- 
lem [24]. Secondly, we can also expand the mode 
functions using an 'initial mode' basis which lead to 
near identical results when compared to the 'instanta- 
neous' method (these were found to be identical when 
there are no intermode couplings, Mmn = 0) [14] 
. We should also include dissipative effects, because 
Im[£i] 7^ for both SPPs and photons. However, 
for SPPs this diminishes the propagation length via 
L = (2Im[A:^])~^, not the production rate. There are 
other issues relating to de-tuning of the resonant fre- 
quencies arising from dissipation, but these are out of 
the scope of the current letter and will be addressed 
elsewhere. 

Finally, given that there are no intermode cou- 
plings (A^mn=0) these SPPs behave as 'pure' modes 



and, for example, have well defined squeezing param- 
eters, which could have applications in quantum in- 
formation. It would also be interesting to investigate 
the vacuum excitation of volume plasmon polaritons 
(VPPs). Due to their longitudinal nature they are not 
excited by light and hence require particle impact at 
the interface, e.g., see [17]. However, vacuum excita- 
tions might be achieved dynamically by firing clusters 
of 'neutral' Argon atoms at a sample of material [25], 
or by placing the sample on a high frequency piezo, 
e.g., see [26]. We hope to comment on this and related 
issues in future work. 

We thank Professor Y. Kido (Ritsumeikan Univer- 
sity), Dr. R. Johansson (RIKEN) and members of the 
Particle Physics Group (Osaka University) for useful 
discussions. 
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